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ABSTRACT 


The  mean  in-line  force  and  the  time  dependent  transverse 
force  acting  on  a  circular  cylinder  undergoing  periodic  trans- 
verse oscillations  in  an  otherwise  steady  flow  was  measured 
for  various  amplitudes  and  frequencies  of  oscillation  at  several 
steady  flow  velocities.   The  experiments  were  carried  out  in  a 
recirculating  water  tunnel  operating  as  a  closed  channel  at 
the  test  section.   The  mean  in-line  force  has  been  expressed 
in  terms  of  a  mean  drag  coefficient  C^i  and  plotted  as  a  func- 
tion of  A/D  and  D/VT .   The  time  dependent  transverse  force  has 
been  expressed  in  terms  of  the  Fourier-averaged  drag  and  inertia 
coefficients  C^  and  Cmi  and  plotted  as  functions  of  the  rela- 
tive amplitude  A/D  and  the  period  parameter  VT/D.   The  results 
have  shown  that  the  mean  in-line  force  is  significantly  larger 
than  that  corresponding  to  steady  flow  about  a  non-oscillating 
cylinder  and  that  energy  may  be  transferred  to  the  oscillating 
cylinder  from  the  fluid  at  or  near  the  vortex-shedding  fre- 
quencies for  the  A/D  values  tested. 
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A  amplitude  of  cylinder  oscillation 

C  ., .  mean  in-line  drag  coefficient 
di  ^ 

C,,  Fourier-averaged  drag  coefficient  (Lift  direction) 

C  ,  Fourier-averaged  inertia  coefficient  (Lift  direction) 

D  diameter  of  test  cylinder 

F  instantaneous  total  force  acting  on  test  cylinder 

T  period  of  oscillation 

t  time 

U  instantaneous  velocity 

U  maximum  velocity  in  a  cycle  (U  =2ttA/T) 

m  2  J  m 

V  mean  velocity  of  the  uniform  flow 

A/D  dimensionless  amplitude  ratio 

D/VT  dimensionless  frequency  parameter 

VT/D  dimensionless  period  parameter 

Re  Reynolds  number  (Re=VD/v) 

6  phase  angle 

v  fluid  kinematic  viscosity 

p  fluid  density 
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I.   INTRODUCTION 

Elastic  structures  of  one  or  more  degrees  of  freedom  can 
extract  energy  from  the  flow  about  them  and  can  develop  cata- 
strophic f low-induced  vibrations.   The  understanding  of  this 
energy-extraction  process  is  of  paramount  importance  if  one 
is  either  to  eliminate  or  minimize  it  or  to  design  the  elastic 
structure  such  that  it  can  withstand  the  oscillations  under 
the  contemplated  environmental  conditions. 

Even  a  superficial  familiarity  with  the  parameters  involved 
in  this  type  of  phenomenon  shows  that  the  investigation  of  the 
flow- induced  vibrations  or  the  flow  about  oscillating  bodies 
is  more  than  an  extension  of  the  past  studies  on  flow  about 
bluff  bodies  at  rest.   In  fact,  the  complexity  of  the  problem 
is  increased  by  an  order  of  magnitude.   In  view  of  this  fact, 
it  is  no  wonder  that  the  past  decade  has  produced  either  num- 
erous experimental  data  with  the  goal  of  obtaining  ad-hoc 
solutions   for  specific  problems,  while  directing  little  if 
any  attention  toward  elucidation  of  the  underlying  mechanics 
and  toward  generalization  of  the  results  for  future  applica- 
tions; or  numerous  analytical  models  most  or  all  of  which  had 
nothing  to  do  with  the  motion  of  the  very  medium  which  pro- 
vided the  necessary  energy  to  set  the  body  in  motion.   In  fact, 
the  fluid  mechanics  of  the  flow-induced  oscillations  became 
so  incidental  to  the  phenomenon  that  the  models  dealt  essenti- 
ally with  black-box-induced  oscillations.   To  be  sure,  these 
models  heavily  relied  on  experimental  data  partly  to  justify 


their  existence  and  partly  to  assign  numerical  values  to 
numerous  variable  constants  imbedded  in  them.   The  ability 
of  these  models  to  scale,  i.e.,  to  permit  extrapolation,  is 
uncertain.   Thus,  in  the  final  analysis  one  is  not  quite 
sure  whether  one  should  use  the  experimental  data  within  the 
range  of  their  application  or  the  curves  fitted  to  them  by 
the  empirical  models.   Suffice  it  to  say  that  the  phenomenon 
is  far  from  understood  even  for  the  idealized  conditions 
encountered  in  the  laboratory  without  the  alarming  consequences 
of  the  real  ocean  environment  where  fauna  and  flora,  ever- 
changing  ocean  currents,  and  temperature  gradients  add  further 
complications  to  an  already  complex  problem.   Surely,  over 
design  is  not  the  answer  but  it  may,  in  the  next  decade  or  so, 
be  the  best  available  tradeoff  with  failure. 

Mathematical  models  of  flow-induced  vibrations  of  bluff 
bodies  and  the  response  of  circular  cylinders  to  vortex  shedding 
have  been  aptly  described  by  Parkinson  [1]  and  Currie,  et  al . 
[2]  and  will  not  be  repeated  here. 

It  appears  that  among  the  various  models  considered  so  far, 
the  "wake-oscillator"  model  of  Hartlen  and  Currie  [3]  attracted 
more  attention  among  the  students  of  vibration  analysis.   We 
will  not  elaborate  here  on  the  attempts  made  by  others  to 
introduce  one  or  more  additional  terms  into  the  model  origi- 
nally proposed  by  Hartlen  and  Currie  for  such  naive  attempts 
produced  only  more  papers  and  unrealistically  defined  coeffici- 
ents than  sound  information  toward  the  understanding  of  the 
flow-induced  oscillations  both  in-line  with  and  transverse  to 
the  ambient  flow. 
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It  suffices  to  conclude  that  despite  the  extensive  atten- 
tion given  to  the  study  of  the  transverse  oscillations  of 
bluff  bodies  in  air  and  water,  there  has  been  very  little  or 
no  attempt  to  measure  the  forces  acting  on  the  bodies  and  to 
separate  the  time-dependent  force  into  suitable  components. 
It  has,  however,  been  firmly  established  that  there  is  a  band 
of  frequencies  in  the  vicinity  of  the  Strouhal  frequency  in 
which  the  vortex  shedding  is  synchronized  or  locked-in  and 
that  both  the  intensity  of  the  vortex  shedding  and  the  span- 
wise  correlation  are  enhanced.   The  synchronization  or  lock-in 
is  manifested  by  a  rapid  decrease  in  inertial  force  and  a 
rapid  increase  in  the  absolute  value  of  the  drag  force.   In 
other  words,  the  total  force  which  is  nearly  in-phase  with  the 
motion  before  synchronization  becomes  nearly  out  of  phase  at 
and  after  synchronization. 

The  purpose  of  the  present  investigation  is  the  determina- 
tion of  these  in-phase  and  out-of-phase  components  of  the  time 
dependent  force  acting  on  cylinders  undergoing  transverse 
oscillations  in  a  uniform  stream  and  the  correlation  of  the 
normalized  force  components  with  appropriate  governing  param- 
eters through  the  use  of  an  equation  proposed  herein.   Evi- 
dently, it  is  through  the  understanding  of  the  forces  acting 
on  oscillating  bodies  that  one  would  eventually  be  able  to 
devise  models  to  predict  the  strumming  or  the  hydroelastic 
behavior  of  rigid  or  elastic  bodies  immersed  in  fluid  streams. 
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II .   EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A.   EQUIPMENT 

1.  NPS  Water  Tunnel 

The  experiments  were  performed  in  a  recirculating  water 
tunnel  (Figure  1)  which  had  a  capacity  of  approximately  500 
gallons.   The  galvanized  test  section,  four  inches  wide,  eight 
inches  high  and  sixteen  inches  long,  was  closed  on  top  with  a 
removable  plexiglass  plate  to  eliminate  the  free  surface 
effects.   A  small  space  adjacent  to  the  side  walls  was  provided 
to  allow  the  passage  of  thin  leaf  arms  connecting  the  cylinder, 
in  the  test  section,  to  the  driving  hardware  above  the  test 
section.   A  low  RPM,  high  capacity,  fourteen-inch-diameter- 
discharge  centrifugal  pump  was  used  to  circulate  the  fluid 
through  the  test  section.   The  velocity  of  the  fluid  was  regu- 
lated by  a  butterfly  valve  arrangement  situated  downstream  of 
the  test  section.   Velocities  of  0.8  to  1.4  feet-per-second 
were  obtained  by  adjusting  the  vanes  of  the  butterfly  valve. 

2.  Motor,  FlywheeT  and  Pivot  Assembly 

The  periodic  motion  in  this  experiment  was  obtained 
by  the  use  of  a  small,  variable  speed,  electric  motor,  a  fly- 
wheel and  pivot  assembly  (Figure  2) .   The  electric  motor, 
mounted  so  as  to  isolate  its  vibrations  from  the  remainder  of 
the  test  apparatus,  transmitted  rotary  motion  through  a  flex- 
ible belt  to  a  flywheel.   A  56  inch  long,  rigid  bar,  pivoted 
at  its  center,  transmitted  the  vertical  component  of  motion 
of  a  bearing,  attached  to  the  flywheel,  to  a  shaft  connected 
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Fig.  2     Motor,  Flywheel   and  Pivot  Assembly 
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to  the  yoke  assembly  and  cylinder.   Horizontal  motion  was 
expended  at  both  ends  of  the  pivot  rod  by  the  free  motion  of 
bearings  along  the  pivot  arm.   The  amplitude,  "A" ,  of  the 
vertical  periodic  motion  was  set  by  adjusting  the  radial  posi- 
tion of  the  bearing  attached  to  the  flywheel.   For  these 
experiments,  the  amplitude  ranged  from  0.25  to  0.75  inches. 
The  frequency  of  oscillation  was  regulated  up  to  four  cycles 
per  second  by  adjustment  of  the  variable  speed  motor.   The 
system  was  designed  and  constructed  to  produce  essentially 
sinusoidal  oscillations,  free  of  secondary  oscillations. 

3.  Test  Specimen  and  Yoke  Assembly 

A  one-inch  diameter  circular  cylinder,  constructed  of 
aluminum  tubing,  was  used  as  the  test  specimen.   The  cylinder 
was  held  in  the  test  section  by  a  yoke  assembly  which  was 
connected  to  the  pivot  arm  by  means  of  a  vertically  constrained 
rod.   The  action  of  the  pivot  arm  caused  the  cylinder  to  oscil- 
late vertically  in  the  test  section,  transverse  to  the  stream 
flow.   The  aluminum  yoke  assembly/force  transducer  (Figure  3) 
was  instrumented  with  the  strain  gages  and  accelerometers 
required  to  monitor  the  forces  and  the  acceleration  felt  by 
the  cylinder. 

4 .  Sensors  and  Sensor  Support 

The  mean  velocity  of  the  fluid  was  derived  from  the 

pressure  differential  of  a  pitot  tube  installed  in  the  test 

section  upstream  of  the  test  cylinder.   One  channel  of  a  two 

channel  recorder  was  used  to  continuously  monitor  the  stream 

velocity  during  experiments.   The  calibration  of  the  pitot 

) 
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Fig.  3  Aluminum  Yoke— Force  Transducer 
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tube  and  transducer  was  accomplished  by  means  of  a  water  bottle 
and  precisely  measured  plates  which  were  used  to  produce  a 
known  water  height  differential. 

The  instantaneous  acceleration  of  the  pivot  arm  was 
sensed  by  an  accelerometer  mounted  on  a  yoke  assembly.   The 
acceleration  signal  was  processed  through  a  filter,  to  remove 
high  frequency  vibrations,  and  was  recorded  on  one  channel  of 
a  two  channel  recorder. 

Instantaneous  forces  were  sensed  by  four  piezo- 
resistive  strain  gages.   A  cantilever  beam  type  arrangement, 
forming  the  top  of  the  yoke  assembly,  as  seen  in  figure  3, 
provided  the  mounting  location  for  the  strain  gages.   The 
forces  sensed  by  this  force  transducer  system  were  recorded 
on  the  two  channel  recorder.   Calibration  of  the  force  trans- 
ducer system  was  accomplished  by  hanging  loads  at  the  mid- 
section of  the  cylinder,  supported  at  the  correct  orientation 
so  as  to  calibrate  for  vertical  forces  for  the  transverse  or 
lift  force  measurements,  and  then  for  the  horizontal,  or  in- 
line force  measurements. 

B.   PROCEDURE 

System  checkout  and  sensor  calibration  preceded  the 
experimental  runs.   Several  test  runs  were  made  to  verify  the 
proper  operation  of  the  test  system,  sensors,  processing  and 
recording  equipment,  and  to  establish  the  test  run  procedures 
to  be  used.   It  was  verified  that  the  mechanical  system  pro- 
duced the  required  periodicity  of  oscillation  for  frequencies 
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below  about  three  cycles  per  second  and  amplitudes  below  one 
inch.  It  was  decided  that  the  system  would  not  be  driven  at 
higher  frequencies  or  larger  amplitudes  during  the  test  runs 
because  of  the  presence  of  excessive  vibrational  effects  be- 
yond these  limits.  It  was  verified  that  the  sensors  and  pro- 
cessing equipment  had  adequate  sensitivity  and  phase  reproduction 
to  allow  measurement  of  the  required  physical  parameters. 

The  first  set  of  experiments  has  as  an  objective,  the 
measurement  of  the  mean  fluid-induced  force  on  the  cylinder,  in 
the  direction  of  the  stream  flow,  (in-line  force)  for  various 
frequencies  and  amplitudes  of  cylinder  oscillation  transverse 
to  the  stream  flow.   It  was  verified  by  experiment  that  there 
was  no  in-line  component  of  inertial  force. 

The  physical  recording  of  the  in-line  force  was  accom- 
plished with  the  Honeywell  recorder.   Pitot  differential  pres- 
sure was  recorded  simultaneously  to  verify  the  magnitude  and 
constancy  of  fluid  velocity  during  the  runs.   Accelerometer 
output  was  recorded  on  the  two  channel  recorder  to  provide 
data  on  oscillation  period  and  phase  angle.   The  flow  velocity 
was  set  at  0.84  feet  per  second,  the  amplitude  at  0.25  inches, 
and  the  data  runs  were  taken  as  the  oscillation  frequency  was 
varied  from  zero  to  about  three  cycles  per  second.   The  charts 
were  annotated  with  run  data  for  later  correlation.   The  pro- 
cedure was  repeated  for  amplitudes  of  0.5  and  0.75  inches  and 
the  entire  series  was  repeated  again  at  a  flow  velocity  of  1.3 
feet  per  second. 
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Subsequently,  the  mean  in-line  force,  the  flow  velocity, 
the  oscillation  period  and  amplitude,  and  the  cylinder  diam- 
eter were  recorded  from  the  charts.   A  mean  drag  coefficient 

C,.  was  manually  calculated  from  this  data  for  correlation 
di  J 

with  the  frequency  parameter  D/VT  and  the  nondimensionalized 
amplitude  A/D. 

The  second  set  of  experiments  had  as  an  objective  the 
measurement  of  the  time  dependent  forces  in  the  direction  of 
cylinder  motion,  i.e.,  transverse  to  the  direction  of  stream 
flow.   Since,  in  this  case,  there  was  an  inertial  force  in  the 
direction  of  interest,  it  was  necessary  to  adjust  the  experi- 
mental procedure  in  order  to  identify  the  instantaneous  value 
of  the  inertial  force.   Separation  of  the  inertial  force  from 
the  total  force  measured  with  the  cylinder  submerged,  was 
required  to  ascertain  the  resistive  force  produced  by  the 
fluid  alone. 

The  physical  recording  of  the  transverse  force  was  accom- 
plished using  one  channel  of  the  two  channel  recorder.   Accel- 
eration was  measured  simultaneously  on  the  other  channel.  The 
water  level  in  the  tunnel  was  alternately  raised  and  lowered 
for  each  run  in  order  to  allow  recording  of  the  total  "Wet" 
force  and  the  inertial  "Dry"  force  at  a  common  frequency  and 
amplitude  of  oscillation.   The  acceleration  trace  provided 
the  reference  for  the  time  correlation  of  the  wet  and  dry 
force  data.   Because  of  the  necessity  to  properly  correlate 
wet  and  dry  phase  information,  extreme  care  was  taken  in  veri- 
fying the  zeros  for  the  acceleration  and  force  recordings. 
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Since  the  time  dependent  force  was  of  interest,  the  mean 
force  and  acceleration  was  used  as  the  respective  trace  zeros. 
This  value  was  found  by  electrically  averaging  the  force  and 
acceleration  signals,  and  setting  the  resultant  mean  value  to 
a  convenient  position  on  the  recording  trace.   The  flow  vel- 
ocity was  initially  set  at  0.84  feet  per  second,  the  amplitude 
at  0.25  inches,  and  runs  were  taken  for  various  oscillation 
frequencies  up  to  about  three  cycles  per  second.   The  procedure 
was  repeated  for  amplitudes  of  0.5  and  0.75  inches  and  then  the 
entire  series  was  repeated  at  a  new  flow  velocity  of  1.33  feet 
per  second.   Figure  4  shows  the  force  and  acceleration  traces 
for  a  typical  experimental  run. 

From  traces  similar  to  figure  4,  the  transverse  force  was 
read  and  recorded  from  both  the  wet  and  dry  run  curves.   Start- 
ing from  the  zero  acceleration  time  (illustrated  on  figure  4) 
the  wet  and  dry  forces  were  recorded  each  0.01  second,  for  at 
least  three  complete  cycles.   This  data  was  then  punched  on  IBM 
computer  cards  for  evaluation. 

The  computer  program  for  calculating  the  transverse  drag 
and  inertia  coefficients,  given  in  Appendix  A,  was  also  used, 
along  with  visual  inspection,  to  determine  the  goodness  of  the 
data  and  the  accuracy  of  the  data  retrieval  procedures.   Com- 
puted values  of  the  mean  wet  and  dry  forces,  if  significantly 
different  from  zero,  indicated  possible  errors  in  the  establish- 
ment of  system  zeros  or  the  presence  of  secondary  undesirable 
effects  from  the  mechanical  system.   After  completion  of  this 
data  assessment,  the  drag  and  inertia  coefficients  were 
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Fig.  4     Force  and  Acceleration  Traces 
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calculated  and  recorded  as  a  function  of  the  period  parameter 
VT/D  for  the  three  values  of  the  nondimensionalized  amplitude 
A/D  investigated. 
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III.   DATA  ANALYSIS 

The  numerical  calculations  as  well  as  measurements  in 
time-dependent  flow  yield  the  resultant  force  as  a  function  of 
time  for  a  given  set  of  numerical  values  of  the  independent 
parameters.   Thus  it  is  not  possible  without  a  suitable  hypo- 
thesis, to  express  the  force  both  as  a  function  of  time  and 
remaining  independent  parameters  as  one  ordinarily  would  in  a 
closed  form  solution.   Such  a  working  hypothesis  is  particu- 
larly necessary  for  the  cable  strumming  problem  since  the 
results  are  to  be  incorporated  into  the  dynamics  of  the  cylin- 
der motion  in  the  form  of  a  forcing  function.   It  should  be 
stated  at  the  outset  that  there  is,  at  present,  no  generally 
accepted  hypothesis  to  decompose  the  time  dependent  force  into 
suitable  components.   As  it  will  be  seen  shortly,  even  the 
existing  hypotheses,  such  as  the  so-called  Morison's  equations, 
are  not  applicable  to  periodic  flows  with  a  non-zero  mean 
velocity. 

Stokes,  in  a  remarkable  paper  on  the  motion  of  pendulums, 
showed  that  the  expression  for  the  force  on  a  sphere  oscilla- 
tion in  an  unlimited  viscous  fluid  consists  of  two  terms,  one 
involving  the  acceleration  of  the  sphere  and  the  other  the 
velocity.   This  analysis  shows  that  the  inertia  coefficient 
is  modified  because  of  viscosity  and  is  augmented  over  the 
theoretical  value  valid  for  irrotational  flow.   The  drag  co- 
efficient associated  with  the  velocity  is  modified  because  of 
acceleration,  and  its  value  is  greater  than  it  would  be  if  the 
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sphere  were  moving  with  constant  velocity.   In  general,  the 
force  experienced  by  a  bluff  body  at  a  given  time  depends  on 
the  entire  history  of  its  acceleration  as  well  as  the  instan- 
taneous velocity  and  acceleration.   Thus,  the  drag  coefficient 
in  unsteady  flow  is  not  equal  to  that  at  the  same  instantaneous 
velocity  in  steady  flow.   Neither  is  the  inertia  coefficient 
equal  to  that  found  for  unseparated  potential  flow.   As  yet  a 
theoretical  analysis  of  the  problem  for  separated  flow  is 
difficult  and  much  of  the  desired  information  must  be  obtained 
both  experimentally  and  numerically.   In  this  respect,  the 
experimental  studies  of  Morison  and  his  co-workers  [4]  on  the 
forces  on  piles  due  to  the  action  of  progressive  waves  have 
shed  considerable  light  on  the  problem.   The  forces  are  divided 
into  two  parts,  one  due  to  the  drag,  as  in  the  case  of  flow 
of  constant  velocity,  and  the  other  due  to  the  acceleration 
or  deceleration  of  the  fluid.   The  concept  necessitates  the 

introduction  of  a  drag  coefficient  C ,  and  an  inertia  coeffici- 

r  d 

ent   C      in    the    expression    for    force.       In   particular    if   F    is    the 
m 

force   per   unit    length    experienced  by    a   cylinder,    then 

2 

F    =    0.5    C,pD  |u|  U      +    C       ptt?-      ~  (1) 

all  m        4        dt 

where  U  and  dU/dt  represent  respectively  the  undisturbed 
velocity  and  the  acceleration  of  the  fluid. 

On  the  basis  of  irrotational  flow  around  the  cylinder, 
C   should  be  equal  to  2  (cylinder  at  rest,  the  fluid  acceler- 
ating; otherwise  C   =1),  and  one  may  suppose  that  the  value 
of  C -.    should  be  identical  with  that  applicable  to  a.  constant 
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velocity.   However,  numerous  experiments  show  that  this  is 

not  the  case  and  that  C-.  and  C   show  considerable  variations 

d      m 

from  those  just  cited  above.  Even  though  no  one  has  suggested 
a  better  alternative,  the  use  of  the  Morison's  equation  gave 
rise  to  a  great  deal  of  discussion  on  what  values  of  the  two 
coefficients  should  be  used.  Furthermore,  the  importance  of 
the  viscosity  effect  has  remained  in  doubt  since  the  experi- 
mental evidence  published  over  the  said  period  has  been  quite 
inconclusive. 

The  drag  and  inertia  coefficients  obtained  from  a  large 
number  of  field  tests,  as  compiled  by  Wiegel  [5],  show  exten- 
sive scatter  whether  they  are  plotted  as  a  function  of  the 
Reynolds  number  or  the  so-called  period  parameter  U  T/D.   The 

reasons  for  the  observed  scatter  of  the  coefficients  C   and 

m 

C.  remained  largely  unknown.   The  scatter  was  attributed  to 
several  reasons  or  combinations  thereof  such  as  the  irregu- 
larity of  the  ocean  waves,  free  surface  effects,  inadequacy 
of  the  average  resistance  coefficients  to  represent  the  actual 
variation  of  the  nonlinear  force,  omission  of  some  other 
important  parameter  which  has  not  been  incorporated  into  the 
analysis,  the  effect  of  ocean  currents  on  separation,  vortex 
formation,  and  hence  on  the  forces  acting  on  the  cylinders, 
etc. 

The  most  systematic  evaluation  of  the  Fourier-averaged 
drag  and  inertia  coefficients  has  been  made  by  Keulegan  and 
Carpenter  [6]  through  measurements  on  submerged  horizontal 
cylinders  and  plates  in  the  node  of  a  standing  wave,  applying 
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theoretically  derived  values  for  velocities  and  accelerations. 
Additional  measurements  have  been  made  by  Sarpkaya  [7]  of  the 
in-line  as  well  as  transverse  forces  acting  on  cylinders  and 
spheres  in  a  sinusoidally  oscillating  fluid  and  it  was  found 
that  the  drag  coefficient  as  well  as  the  inertia  coefficient 
for  a  strictly  sinusoidally  oscillating  fluid  (no  mean  vel- 
ocity) is  a  function  of  U  T/D  and  that  the  effect  of  the 

m 

Reynolds  number  is  rather  secondary  and  certainly  obscured 
by  the  excellent  correlation  of  the  data  with  the  period  param- 
eter U  T/D. 

m  ' 

On  the  basis  of  the  above  discussion,  one  would  assume 
that  Morison's  equation  would  apply  equally  well  to  periodic 

flow  with  a  mean  velocity  where  u  =  V  -  U  cos  6  and  that  C,, 

1  m  dl 

and  C  ,  will  have  constant,  time-invariant,  Fourier  or  least- 
squares  averages.   This,  in  turn,  implies  that  C-,-.  and  C  , 
are  independent  of  the  associated  flow  phenomena.   There  is, 
however,  no  a  priori  assurance  in  the  principles  of  fluid 
mechanics  of  theory  of  models  that  this  is,  in  fact,  the  case. 
Thus  the  effect  of  the  combination  of  a  uniform  current  and 
harmonic  oscillations  on  the  time-average  and  oscillatory 
forces  acting  on  circular  cylinders  will  have  to  be  re-examined 
and  the  limits  of  application  of  the  Morison's  equation  be 
delineated. 

It  is  a  priori  evident  that  both  u  =  -U  cos G  and 
^  m 

u  =  V  -  U  cos  6  yield  the  same  acceleration  du/dt.   Thus,  the 
force  in-phase  with  the  acceleration  in  Morison's  equation 
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remains  unaffected  by  the  presence  of  the  mean  flow.   The 
results  presented  herein  show  that  this  is  not  the  case. 
Furthermore,  the  use  of  the  Morison's  equation  as 


0.5pC, (V-U  cosG) 
d    m 


V-U  cos 

m 


°2  du 

m   4   dt 


(2) 


requires  that  the  time-averaged  drag  force  be  calculated  by 

increasing  the  force  calculated  from  the  steady  flow  by  a 

_  2 

factor  [1  +  0.5(U  /V)   ].    The  results  presented  herein  show 

that  such  an  analysis  appreciably  underestimates  the  measured 
mean  forces.   It  suffices  to  state  that  the  fluid  flow  phen- 
omena for  bluff  bodies  are  significantly  affected  by  the 
combination  of  currents  and  harmonic  oscillations  and  that 
the  results  for  steady  currents  alone  and  oscillations  alone 
cannot  be  combined  to  yield  reliable  estimates  of  forces  due 
to  both  acting  together. 

The  time-dependent  forces  per  unit  length  in  the  present 
study  are  analyzed  according  to  the  following  three-coeffici- 
ent equation: 

2 

F  =  0.5C,pDV2  +  C  p7r§—  ^-  (-U  cos^t)  - 
d        m   4   dt    m   T 


0.5CdDUm  p 


2tt, 
cos— t 


2ir 

cos~-t 


(3) 


which  may  be  written  as, 


„  rF  -,  =  C,  +  C  tt2(U  T/D)  (D/VT)2sin2^t  - 
0. 5PDV2    d    mm''       T 


?         —   2  i    2tt  I    2 n" 
Cd (umT/D)  <D/VT)   cos^-t  cos^-t 


(4) 
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in  which  C   and  C-,  are  given  by  their  Fourier  averages  as 
m      a 

2tt 

C   =  (2U  T/^3D)  /  (FsinG)d0/(U  2D)  (5) 

m      m       J  m 

0 
and 

2tt 
Cd  =  -(3/4)  j   (FCOS  0)dO/(Um2D)  (6) 

0 

Evidently,  C.,  C  ,  and  C  ,  are  functions  of  VT/D  and 
1         a        m       d 

U  T/D  or  A/D.   They  may  depend  also  on  the  Reynolds  number 
which  does  not  explicitly  appear  in  the  above  expression  be- 
cause of  the  assumptions  made  in  the  formulation  of  the  basic 
force  equation. 

In  the  foregoing,  neither  the  coefficient  C.  is  assumed 
to  be  equal  to  the  steady-state  drag  coefficient  for  a  uniform 

flow  at  the  constant  velocity  V,  nor  C   and  C,  are  assumed  to 

J  m      d 

be  identical  to  those  obtained  for  a  strictly  harmonic  oscil- 
lation.  In  fact,  the  results  show  that  C,  =  C,  (steady)  only 

for  U   =0,  and  C,  and  C   are  equal  to  those  obtained  for  the 
m  dm       a 

harmonic  oscillation  only  for  VT/D=0. 

The  equation  proposed  above  is  general  enough  to  be  appli- 
cable to  both  in-line  and  transverse  oscillations.  In  the  case 
of  transverse  oscillations,  however,  the  mean  net  force  in  the 
direction  of  oscillation  is  zero,  i.e.,  C.=0,  and  thus,  one  has 

F  (transverse  force)  =      2         -   2sin^t  - 
0.5pDV2L  ml    m  T 


Cdl  (UmT/D)  2  (D/VT)  2  |cos~t 
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cos— t  (7) 


In  the  foregoing,  the  inertia  and  drag  coefficients  have  been 

denoted  as  C  ,  and  C,_  in  order  to  distinguish  them  from  those 
ml       dl 

corresponding  to  in-line  oscillations.   The  subscript  "1" 
carries  the  meaning  of  "Lift"  or  force  in  the  direction  trans- 
verse to  the  stream. 

Ordinarily,  for  a  perfectly  sinusoidal  oscillation  of  the 
cylinder,  the  coefficients  C  ,  and  C,,  would  be  given  by 
equations  (5)  and  (6).   However,  when  the  oscillations  are 
not  perfectly  harmonic,  it  is  relatively  more  accurate  to  use 
the  velocities  and  the  accelerations  encountered  in  the  exper- 
iments rather  than  assuming  them  to  be  simple  harmonic  motions. 
It  is  with  this  objective  in  mind  that  the  equations  (5)  and 
(6)  were  rewritten  as 

c        =    _2T Q^Fsin    9d9 


'ml 


and 


9    2  2tt      ? 

tt^D    pU         fsinz6de 


2TT 

2  f 

■T  QJ  FcosQdQ 


c,,  =  — u  ^ ,    ;  o) 


2PD^  A2  fcos    e    cosOdo 


0 

as  they  would  be  obtained  from  equation  (3)  in  the  usual 
application  of  the  Fourier  analysis.   Evidently,  had  the  oscil- 
lations been  perfectly  harmonic  the  integrals  appearing  in  the 

denominators  of  the  above  equations  would  have  reduced  to 

2tt 


0 


sin20d9    =  tt  (10) 
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and 

2tt 


/cos2C  |cos  o|d0  = 


The  equations  (8)  and  (9)  together  with  equations  (10)  and 
(11)  would  have  reduced  to  equations  (5)  and  (6). 

Since  in  the  present  study  the  periodic  oscillations  were 
not  perfectly  sinusoidal  partly  by  design  in  order  to  obtain 
greater  generality  and  flexibility  in  the  experimentation  and 
partly  due  to  the  constraints  imposed  in  the  design  of  the 
oscillation  mechanism,  it  became  necessary  to  incorporate  into 
equations  (8)  and  (9)  the  exact  form  of  the  oscillations 
encountered  in  the  experiments.   For  this  purpose  the  dry  force, 
which  is  an  exact  representation  of  the  oscillations,  was  norm- 
alized as 

f(6)  =  Z"^ry, , .,  (12) 

F-dry (maximum) 

and  then  the  equations  (8)  and  (9)  were  rewritten  as 

V 

2T nJ      F-f(6)d9 


C„n    =        -^ UJ *     ^i-,~~  (13) 

J 
m 


'ml  2n2T7 


and 


c    _-_*L_ 

Udl  _^_    2    2 

2pDiT    A 


J*     f(9)di 
0 


2tt 

0/F-f  (6  +  TT/2)de  (14) 


/•f  (6    +    tt/2)  I  f  (0    +    tt/2)    d9 
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It  is  easy  to  show  that  equations  (13)  and  (14)  reduce  to 
equations  (8)  and  (9)  or  to  equations  (5)  or  (6)  for  purely 
harmonic  oscillations  for  which  f(0)  =  sin  8.   The  advantage 
of  the  use  of  the  equations  (13)  and  (14)  is  rather  obvious 
for  all  types  of  periodic  oscillations.   Furthermore,  the 
independent  evaluation  of  the  denominators  of  equations  (13) 
and  (14)  and  their  comparison  with  it  and  8/3  respectively  (as 
they  would  have  been  equal  to   had  the  oscillations  been  har- 
monic) gave  an  indication  of  the  deviation  of  the  observed 
periodic  oscillations  from  a  purely  harmonic  oscillation. 

The  force  acting  on  the  cylinder  in  the  in-line  direction 
due  to  the  oscillations  in  the  transverse  direction  is  expressed 
in  terms  of  a  mean  drag  coefficient,  denoted  by  C  . ,  given  as 

—   _  Force  in  the  in-line  direction  c\z.\ 

di  —  ? 

0 . 5  pLDV 

Evidently,  C , .  ,  C  ,  and  Cjn  depend  on  D/VT,  U  T/D  or  2ttA/D, 
J '   di   ml      dl   ^  m 

and  possibly  on  the  Reynolds  number. 

The  experimental  data  are  analyzed  using  the  computer  pro- 
gram given  in  Appendix  A  written  according  to  the  equations 
(13) ,  (14) ,  and  (15)  and  are  plotted  in  terms  of  A/D  and 
either  D/VT  or  VT/D. 
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IV.   DISCUSSION  OF  RESULTS 

The  results  will  be  discussed  in  two  parts.   The  first 
will  be  the  average  in-line  force  acting  on  the  cylinder 
undergoing  forced  periodic  oscillations  in  the  transverse 
direction.   The  second  will  be  the  time  dependent  transverse 
force . 

Evidently,  the  average  in-line  force  is  coupled  with 
secondary  oscillations  due  to  vortex  shedding.   However, 
such  oscillations  are  rather  small  in  both  steady  and  peri- 
odic flows  and  certainly  not  larger  than  about  seven  percent 
of  the  average  force.   It  is  for  this  reason  that  only  the 
average  of  the  in-line  force  acting  on  the  oscillating  cylin- 
der is  presented  here. 

Figures  5,  6,  and  7  show  the  variation  of  the  normalized 
in-line  force  (C , . )  as  a  function  of  D/VT  for  A/D  =0.25, 
0.50  and  0.75  respectively.   Each  figure  represents  the  data 
obtained  with  two  velocities,  namely,  V=0.84  and  V=1.3.   In 
normalized  form  these  velocities  correspond,  for  the  one 
inch  cylinder  used,  to  the  Reynolds  numbers  Re=VD/v  =  7000, 
and  Re  =  10,833. 

Evidently,  the  in-line  force  increases  with  A/D  since 
the  cylinder,  undergoing  transverse  oscillations,  presents 
a  larger  apparent-projected  area  to  the  mean  flow.   This, 
however,  is  only  part  of  the  explanation.   In  addition,  the 
vortex  growth  and  motion  are  affected  by  the  oscillation  of 
the  cylinder  which  in  turn  affect  the  in-line  and  transverse 
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forces  acting  on  the  cylinder.   This  is  evidenced  by  the  fact 
that  the  in-line  force  for  a  given  A/D  increases  at  first, 
reaches  a  maximum,  and  then  decreases  as  D/VT  increases.   A 
simple  minded  calculation  based  on  the  steady  flow  drag  co- 
efficient for  a  stationary  cylinder  and  the  apparent  projected 
area  for  the  in-line  force  F,  which  may  be  written  as 

C.  •  =   ~ —  =   C,  (1  +  2A/D) 

di     ,  -2         ds 
*5PV  DL 

yields  values  which  are  almost  equal  to  the  maximum  values 
given  in  figures  5,  6,  and  7.   It  should  be  noted,  however, 
that  the  phenomenon  is  far  more  complex,  and  that  such  a 
simple  minded  procedure  should  not  generally  be  used,  even 
though  the  results  are  surprisingly  good. 

For  the  purposes  of  comparison,  the  figures  5,  6,  and  7 
are  combined  in  figure  8  by  drawing  mean  lines  through  the 
data  points.   Figure  8  shows  that  the  in-line  force  coeffici- 
ent reaches  its  maximum  at  D/VT  between  0.18  and  0.20.   Ordin- 
arily, the  Strouhal  number  for  a  stationary  cylinder  would  be 
0.22  for  the  Reynolds  numbers  cited  previously,  and  one  would 
expect  that  the  forces  acting  on  the  cylinder  will  undergo 
dramatic  changes  as  the  vortex  shedding  frequency  given  by 
the  Strouhal  number  coincides  with  the  frequency  of  the  cylin- 
der oscillations.   The  present  results  show  that  such  a  syn- 
chronization takes  place  at  a  frequency  slightly  lower  than 
the  Strouhal  frequency. 

Also  shown  in  figure  8  is  the  normalized  amplitude  of 
the  oscillations  in  the  in-line  force  for  A/D=0.5.   As  noted 
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earlier  the  oscillations  are  quite  negligible  relative  to  the 
mean  force  and  certainly  under  7  per-cent.   It  should  be  noted 
that  the  amplitude  of  the  oscillations,  like  the  average  force 
show  an  almost  sudden  increase  in  the  vicinity  of  D/VT  nearly 
equal  0.19  and  remain  at  that  value  for  larger  values  of  D/VT. 
The  occurrence  of  synchronization  as  well  as  the  increase  of 
the  amplitude  of  oscillations  are  shown  most  dramatically  in 
figure  9.   This  figure  was  obtained  by  setting  the  free  stream 
velocity  at  0.84  feet  per  second  and  the  A/D  ratio  equal  to 
0.5.   Then,  beginning  with  the  case  of  the  non-oscillating 
cylinder,  the  frequency  of  the  oscillations  was  gradually 
increased  up  to  about  four  cycles  per  second  and  the  resulting 
in-line  force  was  continuously  recorded.   The  figure  shows 
that  the  in-line  force  increases  rapidly  but  with  very  little 
oscillations  superimposed  on  it.   As  soon  as  the  state  of 
synchronization  is  reached,  the  amplitude  as  well  as  the  fre- 
quency of  the  force  oscillations  increases. 

From  an  engineering  point  of  view  the  significance  of  the 
magnitude  of  the  in-line  force  is  that  a  cylinder  or  cable 
excited  by  the  flow  to  oscillate  in  the  transverse  direction 
may  be  subjected  to  in-line  forces  several  times  larger  than 
assumed  in  its  design.   Furthermore,  the  deflections  caused 
by  the  in-line  force  on  sufficiently  flexible  cylinders  tend 
to  couple  with  transverse  ocsillations  and  not  only  affect 
the  magnitude  of  the  transverse  oscillations  but  also  the 
path  of  the  cylinder  motion.   Thus  it  is  not  uncommon  to  see 
heat  exchanger  pipes  or  chimneys  exhibit  oscillation  patterns 
in  the  plane  normal  to  their  axis. 
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The  time  dependent  transverse  force  is  described,  as 

noted  earlier,  in  terms  of  a  drag  coefficient  C-,-,  and  an 

inertia  coefficient  C  ,  given  by 

2  2      i 

Transverse  Force=C  , pnD  d  (-U  cos8)  -  C,,pDU   cos  6  cos  0 

ml^—  -rr   m         dip  m       J 

Fiqures  10.  11.  and  12  show  C,,  and  C  ,  as  a  function  of 

*  dl      ml 

VT/D  for  A/D  -  0.25,  0.50  and  0.75  respectively.   These  co- 
efficients were  obtained  without  the  C.  term  in  the  general- 
ized Morison  equation  (see  equation  3)  since  the  mean  of  the 
transverse  force  is  zero. 

It  is  seen  from  these  figures  that  important  variations 

in  C,,  and  C  ,  occur  particularly  in  the  vicinity  of  the 
dl      ml       r 

Strouhal  frequency  (here  VT/D=4 . 5  to  7)  where  the  natural 
eddy-shedding  at  the  Strouhal  frequency  is  both  enhanced  and 
correlated  by  the  oscillations. 

The  inertia  coefficient  or  the  normalized  in-phase  com- 
ponent of  the  transverse  force  undergoes  a  rapid  drop  as  the 
frequency  of  the  oscillation  approaches  the  Strouhal  frequency 
from  both  the  upper  and  lower  limits.   In  other  words,  syn- 
chronization or  lock-in  is  manifested  by  a  rapid  decrease 
in  inertial  force  and  a  rapid  increase  in  the  absolute  value 
of  the  drag  force. 

This  fact,  which  has  not  been  recognized  before,  shows 
that  the  lock-in  phenomenon  is  a  phase  transformer.   The 
total  force  which  is  nearly  in  phase  with  the  motion  before 
synchronization  becomes  nearly  out  of  phase  at  or  after  syn- 
chronization.  It  should  be  noted  in  passing  that  the  success 
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of  many  empirical  models  dealing  with  this  type  of  oscillation 
comes  partly  from  the  manipulation  of  the  phase  angle  near 
synchronization.   It  is  now  apparent  that  the  fluid  force  to 
be  used  in  the  equations  expressing  the  motion  of  a  cylinder 
or  cable  should  be  given  by  the  data  presented  herein.   Such 
data  take  care,  not  only  of  the  variation  of  the  phase  angle, 
but  also  the  amplitude  of  the  transverse  force  as  a  function 
of  the  normalized  frequency. 

Figure  13  depicts  an  example  of  the  occurrence  of  syn- 
chronization as  the  period  of  oscillation  is  gradually  de- 
creased.  The  upper  trace  shows  the  phenomenon  as  frequency 
is  increased  from  the  non-oscillating  case  to  just  beyond 
the  synchronization  frequency.   The  lower  trace  provides  an 
exploded  view  of  the  phenomenon  near  the  synchronization  fre- 
quency.  The  transverse  force  sharply  decreases  as  soon  as 
the  point  of  synchronization  is  reached.   Beyond  that  point, 
the  changes  in  the  transverse  force  with  frequency  are  quite 
small.   The  phenomenon  is  reversible  and  one  would  obtain  a 
figure  similar  to  figure  13  if  one  gradually  decreased  the 
frequency.   The  possibility  of  the  occurrence  of  an  hysteresis 
has  not  been  investigated. 

The  drag  coefficient  C, ,  or  the  normalized  out-of-phase 
component  of  the  total  instantaneous  transverse  force  given 
in  figures  10,  11,  and  12  show  that  C,,  becomes  negative  for 
VT/D  values  between  approximately  4.5  and  7  (i.e.  for  norm- 
alized frequencies  between  0.14  and  0.22).   Outside  this 
range  the  drag  is  positive,  thus  in  the  opposite  direction 
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to  the  motion  of  the  cylinder.   Within  the  range  of  VT/D  values 
cited  above,  the  drag  force  is  in  phase  with  the  direction  of 
motion  of  the  cylinder  and  helps  to  magnify  the  oscillations 
rather  than  damp  them  out.   For  this  reason,  the  region  in 
which  C,,  is  negative  is  sometimes  referred  to  as  the  negative 
damping  region.   The  fact  of  the  matter  is  that  this  is  not 
damping  in  the  proper  use  of  the  word  but  rather  an  energy 
transfer  from  the  fluid  to  the  cylinder  via  the  mechanism  of 
synchronization.   The  values  of  VT/D  at  which  C,,  changes  its 
sign  depend  on  A/D  as  seen  in  figures  10,  11,  and  12  even 
though  the  negative  C, ,  region  roughly  lies  within  the  VT/D 
values  of  4.5  and  7.   The  envelopes  of  the  two  zero  crossings 
determine  the  region  of  self  excited  oscillations.   A  precise 
determination  of  such  an  envelope  and  its  dependence  on  the 
Reynolds  number  will  require  experiments  with  additional  A/D 
values  and  other  Reynolds  numbers.   For  purposes  of  the  present 
investigation,  it  suffices  to  note  that  the  maximum  absolute 
value  of  C-. ,  in  the  synchronization  range  decreases  rapidly 
as  A/D  increases.   Field  studies  have  shown  that  synchroniza- 
tion does  not  occur  for  relative  amplitudes  larger  than  unity. 
The  trend  of  the  present  data  is  in  conformity  with  such 
observations . 

Finally,  an  unexpected  and  previously  unknown  observation 

in  connection  with  the  variation  of  C,n  will  be  described. 

dl 

For  normalized  frequencies  (D/VT)  smaller  than  about  0.14 
the  data  yield  positive  drag  coefficients.   Between  0.14  and 
0.22,  the  synchronization  takes  place  and  the  drag  coefficient 
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is  negative  as  noted  above.   Ordinarily  one  would  have  expected 
that  the  drag  coefficient  will  remain  positive  and  continue 
to  increase  with  increasing  frequencies  beyond  D/VT  =  0.22  and 
eventually  reach  a  value  which  would  be  identical  to  that 
obtained  by  oscillating  the  cylinder  in  a  fluid  otherwise  at 
rest.   However,  an  interesting  phenomenon  takes  place  at  fre- 
quencies between  approximately  0.22  and  0.27,  depending  on 
the  A/D  ratio.   For  example  for  A/D=0.75,  C,,  increases 
rapidly  at  D/VT  =  0.26  and  then  decreases  sharply  to  nearly 
zero.   At  first  it  was  suspected  that  this  might  be  due  to  an 
experimental  error.   However,  the  repeatability  of  the  experi- 
ments and  the  observation  of  the  same  phenomenon  at  other 
relative  amplitudes  and  velocities  have  shown  that  there  is 
indeed  a  dramatic  change  in  C,,  at  D/VT  =  0.26  for  A/D  =  0.75, 
at  D/VT  =0.24  for  A/D  =0.5  and  at  D/VT=  0.31  for  A/D  =  0.25. 
It  is  further  noted  from  the  data  presented  herein  that, 
particularly  for  A/D  =  0.25  and  0.5,  C, ,  becomes  once  more 
negative  in  a  narrower  range  of  D/VT  values  (D/VT  from  0.345 
to  0.45  for  A/D  =  0.25)  and  shows  the  existence  of  a  second 
region  of  synchronization.   The  D/VT  value  at  which  C,, 
acquires  its  second  minimum  value  in  the  case  of  A/D  =  0.25 
is  almost  exactly  twice  that  of  the  first  minimum.   The 
occurrence  of  this  second  region  of  synchronization  at  higher 
frequencies  depends  on  the  A/D  value.   The  fact  that  is 
demonstrated  here  is  that  there  is  not  a  single  region  of 
synchronization  and  that  there  is  at  least  one  and  possibly 
more  regions  of  frequency  in  which  synchronization  can  occur. 
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The  narrowness  of  the  regions  of  frequency  in  the  second 
region  of  synchronization  makes  it  rather  difficult  to  observe 
the  phenomenon.   In  fact  one  may  easily  miss  such  a  region  by 
simply  not  taking  smaller  increments  in  frequency.   It  suf- 
fices to  say  that  a  cylinder  may  be  excited  first  at  fre- 
quencies near  the  Strouhal  frequency  and  then  at  the  multiples 
of  the  Strouhal  frequency.   However,  the  largest  energy  trans- 
fer from  the  fluid  to  the  cylinder  occurs  in  the  first  syn- 
chronization region  near  the  Strouhal  frequency. 
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V.   CONCLUSIONS 

The  experimental  investigation  of  the  transverse  oscilla- 
tions of  a  cylinder  in  a  flow  with  an  ambient  mean  velocity 

V  has  yielded  the  force  coefficients  C , . ,  C  ,  and  C,  and 
J  di    ml      dl 

has  shown  that; 

a.  The  mean  flow  has  significant  effects  on  the  force 
transfer  coefficients  and  that  the  result  of  experiments  with 
harmonic  oscillations  in  a  fluid  otherwise  at  rest  are  not 
applicable  to  the  transverse  oscillations  of  a  cylinder  in  a 
uniform  flow; 

b.  It  is  possible  to  excite  transverse  oscillations  for 
A/D  smaller  than  about  unity.   In  a  region  of  D/VT  in  the 
vicinity  of  the  Strouhal  frequency,  the  out-of-phase  compo- 
nent of  the  total  force  becomes  negative  and  some  energy  is 
transferred  from  the  fluid  to  the  cylinder.   The  rate  at 
which  energy  is  transferred  decreases  with  increasing  rela- 
tive amplitudes.   Furthermore  there  is  at  least  one,  and 
possibly  more,  narrower  bands  of  frequencies  at  which  syn- 
chronization occurs. 

c.  Transverse  oscillations  give  rise  to  an  increai. 
drag  force  in  the  direction  of  the  mean  flow.   This  force 
depends  on  A/D  as  well  as  on  D/VT  and  reaches  a  maximum  at 
about  D/VT  =  0.18.   This  value  corresponds  to  a  normalized 
frequency  slightly  below  the  Strouhal  frequency  and  is  well 
within  the  synchronization  region.   From  a  practical  point  of 
view  this  is  a  matter  of  concern  for  sound  design  of  cables 
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and  other  structures  which  may  be  subjected  to  transverse 
oscillations; 

d.   A  meaningful  dynamic  analysis  of  the  vortex  excited 
transverse  oscillations  may  be  carried  out  only  through  the 
use  of  the  force  transfer  coefficients  presented  herein. 
However,  it  is  necessary  that  additional  data  be  obtained  at 
other  amplitudes  and  Reynolds  numbers.   Furthermore,  additional 
experiments  may  also  have  to  consider  the  roughness  of  the 
oscillating  structure. 
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